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We have carefully analysed the potential of future Cosmic Microwave Background (CMB) and 
Large Scale Structure (LSS) measurements to probe neutrino masses. We perform a Fisher matrix 
analysis on a 9-dimensional cosmological parameter space and find that data from the Planck CMB 
experiment combined with the Sloan Digital Sky Survey (SDSS) can measure a neutrino mass of 
0.12 eV at 95% conf. This is almost at the level of the 0.06 eV mass suggested by current neutrino 
oscillation data. A future galaxy survey with an order of magnitude larger survey volume than the 
SDSS would allow for a neutrino mass determination of 0.03-0.05 eV (95% conf.). 
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I. INTRODUCTION 


The absolute value of neutrino masses are very dif¬ 
ficult to measure experimentally. On the other hand, 
mass differences between neutrino mass eigenstates, 
(toi, m2, m3), can be measured in neutrino oscillation ex¬ 
periments. Observations of atmospheric neutrinos sug- 
rest a squared mass difference of 5m? ~ 3 x 10“^ eV^ 
10, H]. While there are still several viable solutions to 
the solar neutrino problem the so-called large mixing an¬ 
gle solution gives by far the best fit with 5m? ~ 5 x 10“® 

eV^ 1,1. 

In the simplest case where neutrino masses are hi¬ 
erarchical these results suggest that mi ^ 0, m2 ~ 
i^TOsoiar, and 77i^ ^ i^matmospheric- If the hierarchy is 
inverted I, g, 1,1,0, 0 one instead finds m3 ~ 0, 
m2 ~ f^matmospheric, and m\ ~ (^matmospheric■ How¬ 
ever, it is also possible that neutrino masses are de¬ 
generate II |T|, 0 0, 0 0, 0 |I|, 0 0, 0, 
mi ~ m2 ~ m3 ::J> (5matmospheric, in which case oscillation 
experiments are not useful for determining the absolute 
mass scale. 

Experiments which rely on kinematical effects of the 
neutrino mass offer the strongest probe of this overall 
mass scale. Tritium decay measurements have been able 
to put an upper limit on the electron neutrino mass of 
2.2 eV (95% conf.) [0. However, cosmology at present 
yields an even stronger limit which is also based on the 
kinematics of neutrino mass. Neutrinos decouple at a 
temperature of 1-2 MeV in the early universe, shortly be¬ 
fore electron-positron annihilation. Therefore their tem¬ 
perature is lower than the photon temperature by a fac¬ 
tor (4/11)^/^. This again means that the total neutrino 
number density is related to the photon number density 
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( 1 ) 


Massive neutrinos with masses m ^ Tq ~ 2.4 x 10“"* 
eV are non-relativistic at present and therefore con¬ 
tribute to the cosmological matter density U, 0, 0 




92.5 eV’ 


( 2 ) 


calculated for a present day photon temperature Tq = 
2.728K. Here, ^ mj, = mi -I- m2 -I- m3. However, because 
they are so light these neutrinos free stream on a scale of 
roughly k ~ 0.03mev^m^ h Mpc“^ [0 0 0. Below 
this scale neutrino perturbations are completely erased 
and therefore the matter power spectrum is suppressed, 
roughly by AP/P ~ 0. 

This power spectrum suppression allows for a deter¬ 
mination of the neutrino mass from measurements of the 
matter power spectrum on large scales. This matter spec¬ 
trum is related to the galaxy correlation spectrum mea¬ 
sured in large scale structure (LSS) surveys via the bias 
parameter, 6^(fc) = Pg{k)/Pm{k). Such analyses have 
been performed several times before H 0 , most re¬ 
cently using data from the 2dF galaxy survey [ 0 0, 0 . 
These investigations find mass limits of 1.5-3 eV, depend¬ 
ing on assumptions about the cosmological parameter 
space. 

In a seminal paper it was calculated by Eisenstein, 
Hu and Tegmark that future CMB and LSS experiments 
could push the bound on the sum of neutrino masses 
down to about 0.3 eV 0. This calculation was done 
using a Fisher matrix analysis on a 6-dimensional cosmo¬ 
logical parameter space and a fairly crude approximation 
of the MAP CMB data. 

In the present paper we discuss the prospects for neu¬ 
trino mass detection in more detail. First we discuss 
the Fisher matrix technique in some detail and then we 
move on to discuss future data sets from CMB and LSS. 
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It turns out to be possible to describe hypothetical new 
experiments in terms of just a few general parameters. 
Furthermore, the issue of parameter degeneracies is dis¬ 
cussed in some detail. 

As discussed above, the main neutrino physics mass pa¬ 
rameter for CMB and LSS is the sum of all neutrino mass 
eigenstates. However, for hierarchical neutrino masses 
this is almost equivalent to the mass of the heaviest mass 
eigenstate, ^rrii, ~ m 3 . In the remainder of the paper 
we therefore use and ^ interchangeably, except 
where otherwise stated. 


II. FISHER MATRIX ANALYSIS 

Measuring neutrino masses from cosmological data is 
quite involved since for both CMB and LSS the power 
spectra depend on a plethora of different parameters in 
addition to the neutrino mass. Furthermore, since the 
CMB and matter power spectra depend on many differ¬ 
ent parameters one might worry that an analysis which 
is too restricted in parameter space could give spuriously 
strong limits on a given parameter. Therefore, it is desir¬ 
able to study possible parameter degeneracies in a simple 
way before embarking on a full numerical likelihood anal¬ 
ysis. 

It is possible to estimate the precision with which the 
cosmological model parameters can be extracted from a 
given hypothetical data set. The starting point for any 
parameter extraction is the vector of data points, x. This 
can be in the form of the raw data, or in compressed form, 
typically the power spectrum (C; for CMB and P{k) for 
LSS). 

Each data point has contributions from both signal and 


noise, x = ^signal + a^noise- If both signal and noise are 
Gaussian distributed it is possible to build a likelihood 
function from the measured data which has the following 
form 11 ^ 

£(0) oc exp , (3) 

where 0 = (O, iLg, rig, t, ...) is a vector describing 
the given point in model parameter space and C(0) = 
(xx"^) is the data covariance matrix. In the following we 
shall always work with data in the form of a set of power 
spectrum coefficients, Xi, which can be either Ci or P{k). 

If the data points are uncorrelated so that the data 
covariance matrix is diagonal, the likelihood function can 
be reduced to £ oc e~^ where 

^max / \ 2 

2 _ X '' \^i,ohs theoryJ / .\ 

is a x^-statistics and A^max is the number of power spec¬ 
trum data points js^. 

The maximum likelihood is an unbiased estimator, 
which means that 

(0)=0o. (5) 

Here 0o indicates the true parameter vector of the under¬ 
lying cosmological model and ( 0 ) is the average estimate 
of parameters from maximizing the likelihood function. 

The likelihood function should thus peak at 0 ~ 0o, 
and we can expand it to second order around this value. 
The first order derivatives are zero, and the expression is 
thus 


= xLn + 



1 

(x{xkY 


dxk dxk 

WiWj 


{Xk, 


obs 


Xk) 


d'^Xk 

dBidOj 




0 ), 


( 6 ) 


where i,j indicate elements in the parameter vector 0 . 
The second term in the second derivative can be expected 
to be very small because (x^^obs — Xk) is in essence just 
a random measurement error which should average out. 
The remaining term is usually referred to as the Fisher 
information matrix 

a^X^ _ 1 dxk dxk , . 

86 , 86 , a(xfe)2 86 , 86 , ' ^ ’ 

The Fisher matrix is closely related to the precision with 
which the parameters, 0^, can be determined. If all free 
parameters are to be determined from the data alone 
without any priors then it follows from the Cramer-Rao 


inequality that 

a{6,) = (8) 


for an optimal unbiased estimator, such as the maximum 
likelihood Q. 

In order to estimate how degenerate parameter i is 
with another parameter, j, one can calculate how a{6i) 
changes if parameter j is kept fixed instead of free in the 
analysis. Starting from the 2x2 sub-matrix 

(9) 


one then finds 


fixed (^0 



( 10 ) 
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We therefore define the quantity 


_ fixed ^ 

am - 


( 11 ) 


as a measure of the degeneracy between parameters i and 

j- 

In the next two sections we discuss how to calculate the 
contributions to the Fisher matrix from LSS and CMB 
data respectively. 


III. LSS DATA 

At present data from the first very large scale precision 
galaxy surveys are becoming available. The 2dF survey 
has measured about 250,000 galaxies and the Sloan Dig¬ 
ital Sky Survey Q is currently in the process of mea¬ 
suring up to 10® galaxy redshifts. 

With surveys of this scale and precision it will be pos¬ 
sible to obtain quite precise limits on the neutrino mass. 
However, even for large scale surveys the problem is that 
massive neutrinos mainly affect the smaller scales where 
there are potential problems with non-linearity. 

At present analyses are usually carried out using an 
effective cut in fc-space at fc = 0.2/i Mpc“^. This cut is 
placed roughly where the quantity A, defined as 

= J:_4^fc3|4p oc epik), (12) 

is equal to 1 in linear theory. Here, P{k) is the power 
spectrum of fluctuations. In most cases this corre¬ 
sponds reasonably well to the point where 

^LinearI ^Linear! reason for the cut-off is then 

that poorly understood non-linear effects begin to domi¬ 
nate the power spectrum at scales smaller than the cut¬ 
off. 

However, for problems involving the detection of neu¬ 
trino masses in the 0.1 eV range one must be even more 
careful. For instance, for = 0.1 eV the power spec¬ 
trum suppression is only AP/P ~ 0.06 so that any effect 
on the power spectrum at this level which is not well 
understood can impair the mass detection ability. 

From simulations Peacock and Dodds were able 
to derive the following approximate relation between the 
linear and non-linear spectra 

^NL = /nl[Al], (13) 



FIG. 1: The function Q{n) described in Eq. (p2|). 
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3.310(1-b n/3)-°-2^‘‘ 
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Qm — -b (1 -b Hm/2)(1 -b Da/TO) 


(17) 

(18) 

(19) 

( 20 ) 


At the large scales we are interested in one finds that 


/\2 _ A 2 

^Non-Linear ^Linear 
^Linear 


QA 


2 

Linear’ 


( 21 ) 


where 


0 = < 22 , 

We take fim = 0.3 and Da = 0.7 as values for the matter 
and vacuum energy densities respectively. 

In the above equations n is the effective spectral index 
of the linear power spectrum. At very small scales n ^ 
—3 and at scales beyond the horizon at matter-radiation 
equality n ^ 1. Fig. 1 shows Q{n). For the scales we are 
interested in n 0 — 1 so that Q ~ 0.2. 

In order to have a reasonably clean interpretation of 
the data a fc-space cut should be made roughly where 


with 


/nlN = X 


/ 1 -b BPx + [Ax^ \ 

\1 + {[Ax]<^g^{n)/[Vx^m^) ’ 


A = 0.482(1-bn/d)"®-®"*^ 
B = 0.226(1 -bn/3)"^'^^® 


(14) 


(15) 

(16) 


I^Non — Liii0£ir Li 


Linear I 


^Linear 


APik)\m, 

m 


(23) 


For a neutrino mass of 0.1 eV this means that 
O^Linearlfccut - 0-06, Or that ~ 0.3. This 

requirement leads to an estimate of kcut — O.l/i Mpc 
which in fact is not far from the fccut = 0.2/i Mpc“^ often 
used in present analyses. 
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A. fc-dependent bias 

In all present parameter estimation analyses it is as¬ 
sumed that the bias parameter, b^{k) = Pg{k)/Pm{k), is 
independent of the scale, k. 

However, many independent simulations find that bias 
is in fact quite strongly scale dependent ® 0 in the 
non-linear regime. In the linear regime bias is expected to 
be constant, and the asymptotic value boo = limk—,ob{k) 
is reached as a scale of roughly k ~ O.I — 0.2h Mpc“^. 
This means that at scales larger than fccut bias should be 
very close to scale-independent, and that we can therefore 
use a single parameter, &, to describe it. 


Instead of Ves we use Aes = (3I4ff/47r)^/^ as the free 
parameter. As discussed in Ref. [0 the SDSS BRG sur¬ 
vey 0 has an effective volume of roughly {lh~^ Gpc)^, 
corresponding to AeS — 620h~^ Mpc. 

Note that the number of independent Fourier modes 
on a given scale, fc, enclosed within the survey volume is 
proportional to I4ff. Therefore it essentially corresponds 
to the factor {21 -\- 1) for the CMB measurements which 
measures the number of m-modes for a given 1. In that 
sense both 14 ^ and (2^ -|- 1) are a measure of the lack of 
ergodicity in the given data set. 

IV. CMB DATA 


B. Mock LSS surveys 

For purposes of parameter estimation the most impor¬ 
tant parameter in galaxy surveys is the effective volume, 
defined as 


I4ff = 


n{v)P{k) 

1 -I- n{v)P{k) 


(fr. 


(24) 


In the above equation n(r) is the selection function. The 
simple interpretation of V is that it is the volume avail¬ 
able for measuring power at wavenumber k. 

In the following we shall assume that the survey is vol¬ 
ume limited, meaning that the selection function is con¬ 
stant throughout the survey volume. If the survey is flux 
limited the selection function is much more complicated. 
In the region where P{k) > I/h, 14® is independent of k 
and equal to the total survey volume. 

Essentially this means that, with certain restrictions, 
we can use just one free parameter, 14 ®, to describe a 
hypothetical galaxy survey. 

It was shown in Ref. |0 that the contribution to the 
Fisher matrix from such a galaxy survey can be written 
as 


Fij — 2 tt 


(91nP(fc) ainP(fc) 




80 , 


w{k)d\nk, (25) 


where the weight-function is w{k) = 14®/A^ and A = 
211 jk. The upper limit of the integral should be taken 
to be kcut, discussed above. In principle the lower limit, 
knun, should be zero but at large scales the assumption 
that P{k) ^ 1/h breaks down. However, by far the most 
of the weight in the above integral comes from k close 
to the upper limit. Therefore, as long as the k where 
P{k) = 1/h is much smaller than fcmax the error from 
taking 14 ® and fcmin = 0 is quite small. 

It should be noted that the above integral expression 
is quite crude. However, it offers a very simple way of 
estimating parameter estimation errors from galaxy sur¬ 
veys. The error arising from it can be of order a factor 2, 
leading to an error in the estimated (7(0^) of order 2^/^. 


The GMB temperature fluctuations are conveniently 
described in terms of the spherical harmonics power spec¬ 
trum 


Cl = {\al,r^\^), (26) 


where 

A 7^ 

— (0, <^) = X! '^lrnYlm{0, 4>) ■ (27) 

Im 

Since Thomson scattering polarizes light there are ad¬ 
ditional powerspectra coming from the polarization 
anisotropies. The polarization can be divided into a curl- 
free {E) and a curl {B) component, yielding four inde¬ 
pendent power spectra: Ct,i, Ce,i, Cb,i and Cet,i- In the 
present work we neglect the curl-free component since it 
is usually tiny compared with the other contributions. 
However, for some purposes, such as the detection of ten¬ 
sor fluctuation modes, it is essential. Altogether we focus 
on the three spectra Ct,i, Ce,i and Cet,i- 

The initial measurement of CMB anisotropies by 
COBE initiated the development of a new gener¬ 
ation of high-precision CMB experiments. At present 
there are data sets on the temperature anisotropy from 
BOOMERANG [||, MAXIMA 0, DASI 0, VS A 
0 , CBI 0, and several other experiments, as well as 
polarization data from DASI 0]. 

In the near future data from the MAP and Planck 
satellites will become available. This data offers full sky 
coverage combined with very small pixel noise due to long 
integration times. 


A. Foregrounds 

There are several sources of foreground contamination 
in the CMB data. Within the galaxy these are primarily 
synchrotron and free-free emission, as well as dust. Ex- 
tragalactic contamination is mainly due to point sources 
and the SZ-effect. 

Several papers deal with the issue of foregrounds for 
future CMB experiments. In the present work we rely on 
the calculations by Tegmark et al. |0 . Their estimate is 
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that for the temperature anisotropies there is little effect 
from foregrounds at the important frequencies around 
100 GHz. The most likely source of contamination at 
high I is from point sources, but this effect is unlikely to 
be dominant for I < 2500. 

For the polarization anisotropy the problem with fore¬ 
grounds is more severe. For E-type polarization the fore¬ 
ground contamination from point sources is likely to be 
dominant already beyond I ^ 1500. 

In conclusion, no CMB experiment is likely to retrieve 
information on the primary spectrum beyond I ~ 2500 
for temperature and I ^ 1500 for polarization. In the 
present work we assume that the E — T cross-correlation 
is subject to the same foreground contamination as the 
i?-polarization. 

Another important foreground issue is the question 
of secondary CMB anisotropy generation from non¬ 
linearity, i.e. weak leasing and the Rees-Sciama effect 
(and of course the SZ effect discussed above). These ef¬ 
fects have been estimated to be extremely small for the 
range of Fvalues we are interested in here (below 2500), 
but can be dominant on smaller scales. 

Finally, there is a possible effect from inhomogeneous 
reionization. In all present CMB parameter estimation 
analyses reionization is treated as a single parameter 
which is either the optical depth due to reionization, r, 
or the redshift of reionization, z^.. However, the reioniza¬ 
tion must to some extent have been inhomogeneous and 
this could affect the arcminute scale CMB anisotropy. It 
has been estimated for instance by Gruzinov and Hu Q 
that this patchy reionization is unlikely to be important 
for parameter estimation, but whether this is really the 
case is not yet clear. 

For lack of a better description we consider reionization 
in the standard homogeneous picture using r as a free 
parameter, and note that this simplification is unlikely 
to cause any serious problems @ . 


B. Mock CMB experiments 

In view of the discussion above we use the follow¬ 
ing prescription for a mock CMB experiment. We as¬ 
sume it to be cosmic variance (as opposed to foreground) 
limited up to some maximum 1-value. This value can, 
however, be different for temperature and polarization 
detection. Therefore a given hypothetical experiment 
can be described by only two free parameters. It, max 
and Ip,max- For all experiments it will be the case that 

^T,max ^ ^P,max- 

In this picture the MAP data will be well described 
by It, max — 1000 and lp,max = 0, and the Planck data 
by It, max — 2500 and Ip,max = 1500. In some sense 
Planck can therefore be regarded as the “ultimate” CMB 
experiment because is measures all of the power spectrum 
parameter space not dominated by foregrounds. 

The contribution to the Fisher matrix from such a 
CMB experiment is then 


lp,n 




dCi,x „ -l/py n 

■Cov [Ci,x,Ci,y)- 


1=2 X,Y 

^T.max 


de. 


de. 




I —If 


where X,Y = T, A, TE. 

The covariance matrices are given by [pl| 


de. 


Cov{Ci,t,Ci,t) = 
GoyiCi,E,Ci,E) = 
Cov{Ci,te,Ci,te) = 
Cov{Ci,t,Ci,e) = 


(21 + 1 ) 
2 


C} 


l,T 


_ 

(21 + 1) 


(29) 

(30) 


(21 + 1 ) 


[CIt + Ci,tCi,e] (31) 


^2 

^LTE 


(32) 


(21 + 1 ) 

Cov{Ci,t,Ci,te) = (21 -I- (33) 

Cov{Ci,e,Ci,te) = (21 -I- (34) 


It should be noted here that this approximation relies 
on the assumption of 47r sky coverage and no pixel noise 
up to the maximum 1. Even though these assumptions 
are not realised in any real experiment they are suffi¬ 
ciently accurate for estimating the parameter estimation 
accuracy of a given experiment. 


V. NUMERICAL RESULTS 

In order to calculate estimated Icr errors on the vari¬ 
ous cosmological parameters we need to apply the Fisher 
matrix analysis to a specific cosmological model. 

We choose as the reference model the generic ACDM 
model with the following free parameters: Glm, the mat¬ 
ter density, Ha, the vacuum energy density, Hb, the 
baryon density, iJo, the Hubble parameter, rig, the spec¬ 
tral index of the primordial perturbation spectrum, r, the 
optical depth to reionization, Q, the spectrum normaliza¬ 
tion, 6, the bias parameter, and the neutrino mass. 
The reference model has the following parameters: = 

0.3, Ha = 0.7, = 0.02, Hq = 70kms-iMpc"\ 

rig = 1, T = 0, and m^, = 0.07 eV. 

A. CMB data only 

Even though the canonical view of cosmological neu¬ 
trino mass determination is that CMB in itself can 
achieve little because the CMB spectrum is insensitive 
to rrii,, experiments such as Planck will in fact be able to 
probe TTij/ precisely. 
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500 1000 1500 2000 2500 
I 


FIG. 2: The change in Ci due to a neutrino mass of 0.07 eV. 
The shaded band is the cosmic variance error of ACijCi = 
1/V21 + 1. 


As an example of the effect of massive neutrinos on 
the CMB we show the change in temperature and E- 
type polarization spectra for a neutrino mass of 0.07 eV 
in Fig. 2. 

In Fig. 3 we show the expected Icr error bar in a de¬ 
termination of rrii, from a hypothetical future CMB ex¬ 
periment, without invoking any other data whatsoever. 

The hatched region beyond Zp,max = 1500 corresponds 
to the region where the CMB signal is dominated by fore¬ 
grounds, and therefore unlikely to be retrieved in any 
measurement. The thick line corresponds to a{mu) = 
0.07 eV, i.e. the central value suggested by current oscil¬ 
lation analyses provided that neutrino masses are hierar¬ 
chical. 

As can be seen from the figure the Planck satellite in 
itself without any additional information could be able to 
provide a marginal detection of a non-zero rrii,. However, 
an unambiguous result from CMB alone seems unlikely 
and in fact the argument should rather be reversed, as 
was emphasized in Ref. |p^ : Even a neutrino mass as 
small as 0.07 eV will have a significant effect on the power 
spectrum at high I and can bias the estimation of other 
parameters unless accounted for. 


B. The addition of LSS data 

Adding data from LSS surveys markedly improves the 
ability to detect a non-zero itIl,. 

In Figs. 4 and 5 we show how the expected Icr error bar 
is changed for the MAP and Planck experiments when 
LSS data is added. 

In the near future there will be data available from 
the Sloan Digital Sky Survey, for which the Bright Red 
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FIG. 3: Estimated la error bars on rrii, from the Fisher ma¬ 
trix analysis using only GMB data. The thick solid line cor¬ 
responds to rui/ = 0.07 eV, and the hatched region above 
Ip,max = 1500 is the region where the CMB signal is domi¬ 
nated by foregrounds. The blank triangle in the upper left 
corner is the region where Ip,max > It, max. 


Galaxy (BRG) survey in our language corresponds to A ~ 
620h~^ Mpc. In Ref. |2^] the expected la error on TOj, 
was calculated for the case of combining the SDSS BRG 
survey and the MAP satellite data, taking the A:-space 
cut to be at fccut = 0.2h Mpc“^. For this specific case the 
error was estimated to be a{mv) — 0.25 — 0.3 eV. From 
our calculation we find for the same case an expected 
error of 0.23 eV, which can be ascribed to the fact that 
our estimated MAP precision is higher than that used in 
Ref. @. 

For the case of Planck -|- SDSS we estimate the Icr er¬ 
ror to be of the order 0.06 eV. This means that adding 
the SDSS data to the Planck data only improves the de¬ 
tection threshold marginally (from 0.07 eV to 0.06 eV). 

In order to achieve a 2cr detection of a neutrino mass 
of 0.06 eV, which is the most favoured value using the 
present data, it is necessary to combine the Planck data 
with data from a galaxy survey with A 3h ^ Gpc (if 
kcut = O.l/i Mpc“^). Such a survey would contain a large 
fraction of the present Hubble volume and is probably not 
feasible. If, on the other hand, numerical simulations 
could provide the ability to disentangle the non-linear 
effects beyond k ^ O.lh Mpc~^ a much smaller survey 
would suffice. For instance if fccut = 0-5h Mpc“^ the 
survey volume needed for a 2cr detection would be “only” 
1500 h-^ Mpc. 




















7 



500 1000150020002500 
X {hr' Mpc) 


500 1000150020002500 
X (hr' Mpc) 


1.0 0.5 0.0 - 0.5 - 1.0 - 1.5 

log C7(m„) 

FIG. 4: The estimated la error bar on using CMB data 
from MAP combined with a hypothetical LSS survey. The 
thick solid line corresponds to niv = 0.07 eV. 
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FIG. 5: The estimated Icr error bar on m^, using CMB data 
from Planck combined with a hypothetical LSS survey. The 
thick solid line corresponds to rrh^ = 0.07 eV. 


C. Other probes 

Finally, we note that all of these calculations were done 
for the case where there is no prior information on any 
cosmological parameters. If the neutrino mass were the 
only free parameter, SDSS+Planck could achieve a Icr ac¬ 
curacy of 0.003 eV, corresponding to a many-sigma de¬ 
tection of a neutrino with mass of order 0.06 eV. This 
is of course wishful thinking in the sense that no other 
experiment is likely to yield such information on the cos¬ 
mological parameters. 

In Fig. 6 we show a matrix of the parameter for the 
case of Planck-bSDSS data. We note that the main pa¬ 
rameter degeneracies in this case are with Hq and fit,, not 
with flm and b as with the present data. The estimated 
Icr errors on these two parameters from the Fisher ma¬ 
trix analysis are AHq/Ho = 0.018 and AOb/flt, = 0.025. 
If it were possible to measure these two parameters more 
precisely by other means the expected precision of the 
neutrino mass determination could be improved by al¬ 
most a factor 2, to a level of 0.03 eV for SDSS-I-Planck. 


VI. DISCUSSION 

We have discussed in detail the prospects for detect¬ 
ing neutrino masses with CMB and large scale struc¬ 
ture observations, extending the pioneering calculation 
in Ref. Q. 


From the Sloan Digital Sky Survey and the Planck 
data we estimated that it would be possible to obtain a 
Icr error bar on the neutrino mass of roughly 0.06 eV. It 
should be noted that this estimate is probably fairly safe 
in the sense that it requires only robust assumptions. 
For instance it requires only LSS data at scales larger 
than k ~ O.l/i Mpc“^ where effects of non-linearity are 
negligible. 

This data would allow for a marginal detection of the 
neutrino mass if the neutrino mass is 0.06 eV, the value 
favoured by current oscillation data. Maltoni et al. 
find that the most massive mass eigenstate should be at 
m ~ 0.04 — 0.07 eV (assuming the LMA Solar solution). 
Therefore it seems unlikely that neutrino masses will be 
detected by SDSS-I-Planck data, unless additional infor¬ 
mation can be retrieved. Perhaps the most promising 
prospect for this is that detailed numerical simulations 
will allow for a disentanglement of the non-linear effects 
i LSS data at high k. If SDSS data up to fc ~ 0.5h Mpc“^ 
could be used the error would diminish to 0.04 eV. 

However, in order to achieve an unambiguous detec¬ 
tion of a neutrino mass of 0.06 eV a much larger survey 
than the SDSS is needed. Whether such a survey be¬ 
comes feasible in the future remains to be seen, but with 
improvement in detector technology and the easier access 
to 8 meter class telescopes it may not be impossible. 

It should also be noted that in the case of an inverted 
mass hierarchy there will be two mass eigenstates of al¬ 
most degenerate mass, so that the atmospheric neutrino 
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FIG. 6: Values of the parameter Vij, defined in Eq. ([li|), cal¬ 
culated for the specific case of Planck-|-SDSS data. 


data indicates a sum of order 0.12 eV. Such a scenario 
would be detectable at the 2cr level with Planck-|-SDSS. 

All of the mass constraints discussed here should be 
compared with the expected precision of other, planned 
experiments. The KATRIN tritium endpoint experiment 
|n3 is planned to measure the effective electron neutrino 

mass, i to an accuracy of 0.35 

eV (95% conf.). This is better than the 95% accuracy 


expected from MAP-l-SDSS (0.45 eV at 95 % conf.), but 
not competitive with the 0.12 eV (95% conf.) expected 
from SDSS-|-Planck. However, it is possible that tritium 
decay experiments can be pushed to the limit of ^ 
0.06 eV, which would indicate hierarchical masses. 

If neutrinos are Majorana particles processes such as 
neutrinoless double beta decay will become possible. The 
detection of 2/30 j^ decays probes the mass combination 


Ulee 


j 


(35) 


and has led to a current upper bound of 0.27 eV 
Future experiments, such as GENIUS Q, could take 
the accuracy to O.Ol eV, much better than the expected 
accuracy from cosmology. However, this requires that 
neutrinos are Majorana particles and that there is no 
cancellation of terms in the expression for mee- 

In conclusion, using LSS and CMB data which will 
become available within the next 8-10 years and reason¬ 
able assumptions it should become possible to achieve a 
95% conhdence limit on the mass of the heaviest neu¬ 
trino mass eigenstate of 0.12 eV. This is almost at the 
level expected if neutrino masses are hierarchical, and 
with a future generation of LSS survey it seems feasi¬ 
ble to detect neutrino masses as low as 0.03-0.06 eV at 
95% confidence. The timescale for such future surveys is 
comparable to, or shorter than, the timescale for building 
tritium endpoint experiments beyond KATRIN. There¬ 
fore there is good reason to be optimistic regarding the 
future potential of cosmology to probe neutrino masses. 
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